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We analyze an asymptotically exact solution for the transition temperature of p-wave supercon- 
ductivity near ferromagnetic criticality on the basis of three-dimensional electron systems in which 
scattering processes are dominated by exchange interactions with small momentum transfers. Taking 
into account all Feynman diagrams in the gap equation, we show that vertex corrections neglected 
in the conventional Eliashberg formalism enhance the dynamical retardation effect of the pairing 
interaction, and increase the superconducting transition temperature significantly, though they just 
give subleading corrections to properties of the normal state. 



The recent discovery of ferromagnetic superconductors such as UGe2, URhGe, and ZrZn2 has brought about 
renewed interest in superconductivity (SC) in the vicinity of ferromagnetic criticality (FC) [1-3]. The proximity to 
FC implies a possible spin- fluctuation- mediated Cooper pairing between parallel spins [4-11]. For more than 20 years, 
the transition temperatures T c of paramagnon- mediated SC have been evaluated by solving the Eliashberg equation in 
which vertex corrections are neglected [4,5], though it has been recognized by researchers that the vertex corrections 
are important in the paramagnon theory [12]. A few studies based on finite-order perturbative calculations have 
revealed that a certain type of vertex correction increases T c [13]. Nevertheless, it is highly desirable to elucidate 
the role of vertex corrections played in spin-fluctuation-mediated SC by non-perturbative approaches. To obtain a 
better physical insight into this issue, in the present paper, we analyze an electronic system with a pure repulsive 
interaction which is exactly proved to exhibit ferromagnetic critical behaviors and also the transition to spin-triplet 
SC. The superconducting transition temperature T c is determined by including full-order vertex corrections, and 
compared with that obtained by neglecting vertex corrections. We will show that although vertex corrections to the 
single-particle normal self-energy are negligible near FC, T c is substantially increased by vertex corrections to the gap 
equation. 

The model considered in the present paper is a three-dimensional (3D) electronic system in which the ferromagnetic 
spin exchange interaction is dominated by forward scattering with small momentum transfers. The model Hamiltonian 
is of the form, 



where S(q) = ^2k c ak+q eTa l 3C k ^ = UQ{q c — \q±_\) with 6(2;) = 1 for x > and for x < 0, and U is positive. 

In this system, momentum transfers in the direction tangent to the Fermi surface q± are restricted to a small range 
\q±\ < q c , though momentum transfers in the direction normal to the Fermi surface are not. The system can be 
regarded as a low-energy effective theory for electron systems near FC. 

Since, in 3D itinerant electron systems close to FC, the Gaussian fixed point is stable [14,15], a RPA-like treatment 
of the ferromagnetic spin fluctuation is valid for our system. Moreover, because of the reason explained below, the 
spin correlators of model (1) are almost exactly given by RPA-like expressions with unrenormalized parameters. It 
is noted that model (1) is a variant of the fermion systems with strong forward scattering that has been studied 
extensively by several authors [16-18]. In these systems, the expansion with respect to a small parameter q c /(2kp) 
enables us to treat electron-electron interaction in all orders in terms of the coupling strength U. Since the low-energy 
properties of the systems are mainly determined by scattering processes with small momentum transfers, the results 
obtained up to leading order in q c /{2kp) are asymptotically exact at sufficiently low temperatures. To implement this 
analysis, we exploit the remarkable theorem proved by Metzner et al. and Kopietz et al.; for leading order in q c j (2kp), 
fermionic loop diagrams with more than two insertions cancel each other in strong forward scattering fermion systems 
that do not involve spin flip processes [16,17]. The proof given by these authors for this loop cancellation theorem 
can be straightforwardly generalized to the case of spin exchange interactions in the paramagnetic state; i.e. the loop 
cancellation theorem holds, even when spin flip occurs at vertices. However, we should keep in mind the following 
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limitations of the theorem. In the ferromagnetic state, the cancellation among loop diagrams involving spin flip is not 
complete, as is easily checked by considering a loop with three legs involving spin exchange processes. Thus in the 
following we mainly concentrate on the paramagnetic state. A more serious restriction of our argument is that the 
theorem is not applicable to the superconducting state below T c , because loop diagrams involving the condensation 
of Cooper pairs do not cancel [19]. Despite these limitations, the theorem is still useful in analyzing the transition 
temperature T c , which is determined by the linearized gap equation as shown below. 

The loop cancellation theorem ensures that the effective interaction between electrons with parallel spins for model 
(1) is the RPA-like form, 

n , , -V 2 (g)xo(q^ n ) _ -3U 

D ^ n) =l-V( q)xo ( q ,u Jn )~ t+{ ^ r + ^- (2) 

Here X o(q,u n ) = -TJ2 km G° k {e m )G° k+q (e m +Lo n ) with G° k (e m ) = l/(ie m - E k ), and t = 3[1/(IW(0)) - 1], JV(0) - 
kp/(ir 2 Ep). v is the Fermi velocity. Also, the effective interaction which involves spin exchange D q ex \u>) is given 
by D q ex \u>) = 2D q {uj) because of the spin rotational symmetry. In a similar manner, we can express exact spin 
correlation functions as RPA forms with unrenormalized parameters. Thus the transition to a ferromagnetic state is 
completely determined by the Stoner condition. In the following, we apply our analysis only to the paramagnetic case 
LW(0)<1. 

In general, RPA expressions with unrenormalized parameters are exact in the case where interactions between 
fluctuations of order parameters are infinitely long-ranged, and thus all self-energy corrections vanish [20,21]. We 
would like to emphasize that, in our system, the absence of self-energy corrections in cq. (2) is not due to this reason, 
but caused by loop cancellations for finite range interactions, and that the single-particle self-energy does not vanish 
by itself in model (1). 
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FIG. 1. (a) Feynman diagrams of three-point vertices. A^ 1 ' is the normal part. is the anomalous part, (b) Diagram of 
the normal self-energy. The wavy line indicates the effective interaction, (c) Diagrams of the anomalous self-energy. 

For the exact analysis of the superconducting transition, we calculate three-point vertices including both normal 
and anomalous parts, whose diagrammatic expressions are shown in Fig. 1(a). As shown by Metzner et al., in systems 
with strong forward scattering, the irreducible current vertex Afc j / c+9 (e,e + w) is related to the irreducible density 
vertex Afe ! fe+ 9 (e, e + w) as \k t k +g (e, e + ui) « v k+q ^ 2 ^k,k+ q {£, £ + w), because of the velocity conservation [16]. Then, 
we obtain the asymptotic Ward identity for the normal part of the density three-point vertex, 

AgU.,.^)- ^"^'-^ . (3) 
k ' k+q iu-q-v k+q/2 

where Gfe(e) is the exact single-particle Green's function. The anomalous part of the density three-point vertex also 
satisfies the generalized asymptotic Ward identity, 

Ai 2 Ue, = - Afc( : £)+Afc+g(£+ " ) - (4) 

' +q ioj-q-v k+q/2 

Here Afc(e) is the anomalous self-energy, which emerges as a result of the condensation of Cooper pairs. 

Using the Ward identities (3) and (4), and eq. (2), we can write down the exact normal self-energy Efc(e) expressed 
diagrammatically in Fig. I (b), 
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^ H e - e ) - ( fc - fc ) ' W(fe+fc')/2 

fei{e>-e)-{k'-k)-v {k+k , )/2 > 1 j 

where the prefactor <y = 3 stems from the SU(2) spin rotational symmetry, and the exact single-particle Green's 
function is Gfc(e) = [Zk{e)ie — Ek}^ 1 with Zk(e) = 1 — Efc(e)/(«e). The anomalous self-energy (e) at T — T c can be 
linearized in terms of the gap function Ak(e) / Zk(e). The exact linearized gap equation expressed diagrammatically 
in Fig. 1(c) is then written as, 

Ak(e) = G-k'(-e')Dk-k'(e-e')A { k%(e,e') 

k',s' 

+ £ A fe < (s')Gk' (e')G-k' (-e')Dk-k' (e - e 1 )^, (e, e') 

k',e' 

= 2 ^ A fe , (e')Gk' (e')G-k> {-e')D k - k > (e - e')A k % (e, e') 

fc'.e' 

+Afe(ej > tj- : 777 tt 

P^, H £ - £ ) - ( k - k ) ■ V(k+k')/2 

_ y> A fc , (e')£> fc - fc , (e - eQG fc > (e')G fc , (-gQG^ 1 (e) 
^ i(e' - e) - (fe' - fe) ■ v {k+ k>)/2 



fc',e' 



Note that the above equations are exact in the sense that full-order Feynman diagrams are taken into account, going 
beyond the conventional Eliashbcrg equations. In eqs. (6) and (7), the effects of superconducting fluctuations, which 
generally decrease T c and give rise to pseudogap phenomena [22,21], do not exist, because the loop cancellation 
suppresses their contributions. Therefore, T c is exactly determined by competition between the pairing interaction in 
(7) and the pair-breaking effect caused by the normal self-energy (6). Since the gap equation should be solved in the 
vicinity of the Fermi surface at low temperatures, G^ 1 (e) in the numerator of the last term of eq. (7) suppresses this 
term by a factor ~ T/Ep compared to the other terms. Thus we drop it in the following. 

To carry out a more explicit analysis, we consider an isotropic system with a spherical Fermi surface, and linearize 
the dispersion relation as E k = vk. Then, we can easily show that the constant C in eq. (6) is equal to zero. Also, 
in this case, (A;' — fe) • V(k+k')/2 ~ (k' r — k r )v + 2kp v sin 2 (#/2), where k r is the component of momentum normal to 
the Fermi surface, and 9 is the angle between k and k. This approximation is adopted in the denominators of the 
vertex functions in eqs. (6) and (7). Note that we retain the term quadratic in the transverse momentum transfers, 
sin 2 (0/2), because, as was noted by Metzner et al. [16], this term is necessary to eliminate an artificial singularity in 
the vertex functions which destroys the Fermi liquid state. In the vicinity of FC, since low-energy fluctuations are 
dominant, we can neglect the momentum dependence of the self-energy, Zk(e) — ► Z(e). We also expand A& into the 
Z-th angular momentum components A;. Then, integrating over E k >, we recast eqs. (6) and (7) into 



[1 - Z(e n )]e n = -gT V / dxxD (x, s n - e m ) 

m J ° 

i(S(e n ) - S(e m ))(sign£ m - sign(e m - £„)) 



signe m + 



\Z(e m )e m \ + \e m - e n \ + 2ik F vx 2 sign(e m - e n ) _ 

A/(e m ) 



(8) 



[1 - (2Z(e„))" 1 ]A / (e„) =TV f dxxD^x, e n - e m )- 

m Jo 

(9) 



1 + 



\Z(e m )e m \ 

z(S(e„) - S(e m ))sign(e„ - s m ) 



Z(e m )e m \ + \e m - e n \ + 2ik F vx 2 sign(s m - e„) 



~ . ,. 3tt NityUPAl - 2x 2 ) 

D l (x,e-e') = - — - ' „ , (10) 

2 S (e-e>)(S(e-e>) + ^)' 
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where S 2 (uj) = t + x 2 + 3it\uj\ / (AvUfx), Pi{x) is the Legendre polynomial, and x c is a cutoff which restricts the integral 
over the transverse momentum within the range q^ < q c . 

We see from eq. (8) that the vertex corrections to the normal self-energy contribute only in low-energy regions, 
because of the factor, signe r „ — sign(£ m — e„), yielding only subleading corrections compared to the first term of the 
right-hand side of cq. (8). Thus, the low temperature behaviors of the normal self-energy are dominated by a one-loop 
diagram without vertex corrections. To check this more explicitly, we analytically continue to real frequencies, and 
obtain the low temperature behaviors of the normal self-energy. For T < T* = t 3/ ' 2 £ , i?/(1.57r 2 ), the imaginary part of 
the self-energy shows the Fermi liquid (FL) behavior, 

ImSjf (0) = - 27 3 fc 2 ^3 7 r 2 )2 + 0{{T/E F f log(t)A). (11) 

Here the second term is the subleading contribution from the vertex corrections. Accordingly, the mass renormalization 
factor is a finite constant, 

For T > T* , strong spin fluctuations give rise to non- Fermi liquid (NFL) behaviors, 

Qk 3 JJ T 

ImS«(0) = log(T/T*) + 0{{T/E F )% (13) 



A = ^log(Hl^) + 0((T/^). (14) 

The leading temperature dependences of the single-particle properties coincide with those obtained by lowest order 
(self-consistent or non-self-consistent) calculations [23,24] and renormalization group calculations [14,15] up to constant 
factors. In the NFL regime, since the enhanced single-particle damping of the same order as the thermal excitation 
energy T breaks the quasi-particle picture, a perturbative expansion up to finite order in terms of the interaction 
strength is insufficient for the description of the low-energy properties. However, the above exact results show that 
contributions from higher-order diagrams including the vertex corrections (the second terms of the right-hand side of 
eqs. (13) and (14)) are negligible for model (1). This observation implies that the classical theory for ferromagnetic 
fluctuation developed many years ago gives a good approximation for our system [23,24]. 

In contrast to the normal self-energy, the anomalous self-energy enjoys significant contributions from the vertex 
corrections, because, as seen from eq. (9), they are not restricted to low-energy regions. These vertex corrections 
enhance the dynamical retardation effect of the pairing interaction. To examine to what extent this effect affects the 
transition temperature T c , we solve the gap equation (9) numerically. For this purpose, as is usually done, we recast 
cq. (9) into the eigenvalue problem, p(T)A;(e„) = ^ m K mn /S.i{e m ). T c is determined by the temperature at which 
p(T c ) = 1. We consider only the p-w&vc pairing, because it gives the highest T c of our model. The superconducting 
T c calculated as a function of the parameter t, which is the measure of ferromagnetic criticality, is shown in Fig. 2 
for a particular set of parameters. Since our analysis is exact up to leading order in q c /(2kp) = x c , the errors of the 
computed T c are of the order x c T c . For comparison, we also show the transition temperatures calculated by neglecting 
the vertex corrections, T c o, in the same figure. The asymptotically exact T c is always higher than T c o, even if we take 
into account the errors of the order x c T c . The maximum ratio of T c to T c o is approximately ~ 3. It is noted that in 
the vicinity of the FC (t ~ 0), the transition temperature decreases as t approaches zero. This tendency is similar 
to the results obtained from the conventional Eliashberg equations [8] . According to Roussev and Millis [8] , T c does 
not vanish at t = within the Eliashberg formalism. Although it is rather difficult to compute T c at t — with 
numerically reliable accuracy from eqs. (8) and (9), we speculate that it is non- vanishing, since the vertrex corrections 
always increases T c . The above results indicate that model (1) is a rare example of a purely repulsive electron system 
that has been exactly proved to undergo the superconducting transition [20] . 
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FIG. 2. T C /E F versus t for x c = 0.3 (•) and x c = 0.2 (x). T c0 /E F versus f for a; c = 0.3 (o) and :r c = 0.2 (□). The broken 
line indicates the crossover temperature T* which distinguishes the NFL state (high temperature region) and the FL state (low 
temperature region). The solid and dashed-broken lines indicate T c x 0.1/Ef and T c q x 0.1/Ef for x c = 0.3 in the Ising case, 
respectively. 



Next we consider the case with the Ising-like anisotropy of the exchange interaction. To simplify the calculation, we 
assume that the Ising anisotropy of spin degrees of freedom is so strong that the interaction part of the Hamiltonian 
(1) is replaced with — J2 q V(q)S z (q)S z (— q). In this case, the asymptotically exact T c is obtained from eqs. (8) and 
(9) with only one modification that the factor g in (8) is replaced with g = 1. The transition temperatures computed 
for this Ising case are shown in Fig. 2 (solid line). The transition temperature without vertex corrections, T c0 , are 
also shown in the figure (dashed-broken line). T c takes a maximum value at the ferromagnetic critical point t = 
as in the case without vertex corrections [6,8,9]. This is due to the absence of the transverse spin fluctuation which 
causes the pair-breaking effect, and suppresses T c for small t in the SU(2) symmetric case. We would like to note that 
in the case with the strong Ising anisotropy, the loop cancellation theorem holds even in the ferromagnetic metallic 
state. We can calculate the asymptotically exact T c in the ferromagnetic state by using eqs. (6) and (7) with g = 1, 
which may be relevant to ferromagnetic superconductors such as UGe2, URhGc, and ZrZn2. 

In conclusion, we have obtained the asymptotically exact solution for the transition temperature of p-wave SC near 
FC. We have shown that vertex corrections neglected in the conventional Eliashberg equation enhance the dynamical 
retardation effect of the pairing interaction, increasing the transition temperature significantly. 
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